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Abstract — This paper addresses the decoupled control of the 
nonlinear, multiinput-multioutput, and highly coupled space sta- 
tion furnace facility (SSFF) thermal control system. Sliding mode 
control theory, a subset of variable-structure control theory, is 
employed to increase the performance, robustness, and reliability 
of the SSFF’s currently designed control system. This paper 
presents the nonlinear thermal control system description and 
develops the sliding mode controllers that cause the interconnect- 
edw subsystems to operate in their local sliding modes, resulting 
in control system invariance to plant uncertainties and external 
and interaction disturbances. The desired decoupled flow-rate 
tracking is achieved by optimization of the local linear sliding 
mode equations. The controllers are implemented digitally and 
extensive simulation results are presented to show the flow- 
rate tracking robustness and invariance to plant uncertainties, 
nonlinearities, external disturbances, and variations of the system 
pressure supplied to the controlled subsystems. 

Index Terms — Decoupled control, nonlinear control, robust 
control, sliding mode control, thermal control. 

I. Introduction 

T HE space station furnace facility (SSFF) provides the 
necessary core systems to operate various material pro- 
cessing furnaces in the microgravity environment of the In- 
ternational Space Station Alpha. The thermal control system 
(TCS) is defined as one of the core systems, and its function 
is to collect excess heat from furnaces and to provide precise 
cold temperature control of components and of certain furnace 
zones. The TCS utilizes single phase water as its cooling 
medium and consists of piping, heat exchangers, coldplates, 
cooling jackets, pressure, flow, and temperature sensors, and a 
centrifugal pump. The TCS schematic is illustrated in Fig. 1. 

The objective of the TCS is to simultaneously control 
the thermal environments of several furnaces and related 
subsystems plumbed in parallel. Specifically, the high accu- 
racy, robust, and decoupled tracking of time-varying flow-rate 
profiles in paths 1, 2, and 3 is desired. The flow-rate through 
path 4 is passively controlled to a constant by fixing its valve 
position and by maintaining the system pressure drop constant. 

The current SSFF control system consists of four linear 
proportional integral derivative (PID) controllers [5]. Three 
of these track flow profiles in paths 1, 2 and 3, respectively. 
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and one controls the system pressure by modulated the pump’s 
speed based on sensed system pressure. System pressure con- 
trol serves to decouple the flow dynamics between respective 
paths, which aids in flow tracking performance. 

The SSFF will remain aboard the Space Station for at least 
1 5 years. As a result of this long duration mission, the loss or 
degradation of pressure control is a likely event. This event 
would cripple flow tracking performance and would result in 
a total SSFF shut down. Reliability of such systems can be 
improved and performance maintained by the use of powerful 
nonlinear control methods. The most competitive methods 
are feedback linearization and sliding mode control. Both 
methods can provide desired linear performance to nonlinear 
systems; however, feedback linearization is more sensitive 
to uncertainties and disturbances and is more complicated 
to implement than sliding mode control. This article shows 
how sliding mode control will improve reliability as well as 
performance, without adding system complexity. 

According to Itkis [1], Utkin [2], and DeCarlo [3], slid- 
ing mode control is unique from other linear and nonlinear 
control methods in its ability to achieve accurate, robust, 
decoupled tracking for a class of nonlinear time-varying multi- 
input-multioutput systems in the presence of disturbances and 
parameter variations. 

Section II formulates the TCS problem including a descrip- 
tion of the nonlinear multiinput-multioutput system mathe- 
matical model with partially known parameters. Section III 
develops the sliding mode control theory in the context of 
the TCS problem. Section IV discusses the sliding mode 
realization. In Section V, the simulation results are presented 
to demonstrate sliding mode's increased tracking performance, 
robustness, and reliability compared with the linear TCS 
controllers. 

II. TCS Problem Formulation 

This section formulates the TCS control problem upon 
which the sliding mode controllers are designed to solve. The 
problem formulation begins with the derivation of the TCS 
mathematical model. 

A. TCS Mathematical Model 

The mathematical model of the TCS is derived as a means 
to simulate the dynamic response of the flows, pressures, and 
actuators within the system. The model is divided into hy- 
draulic, actuator, and pump submodels, where each submodel 
parameter is defined in the Appendix. 
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Fig. 1. TCS schematic. 


1 ) Hydraulic Submodel : Hydraulic systems such as the 
TCS usually consists of liquid-filled tanks connected by 
pipes having orifices, valves, or other flow-restricting devices. 
The hydraulic mathematical model of the TCS is derived in 
terms of three basic elements: resistance, compliance (fluid 
capacitance), and inertance (fluid inertia) [4]. Shown in Fig. 2, 
these elements are combined to form an electrically equivalent 
circuit for the TCS model. From this circuit, the flow and 
pressure equations are written using KirehofPs voltage law. 
The resulting equations are shown as 

t ’=h( ws - , t w ) 

W S = ^-(P S -P- R S W S ) 

Wi = — — [P — (R v . *f Rp.)Wi], Vi = 1,4 (1) 

l p, 

where P = Pi - P 2 is the system pressure drop, Ws is the 

pump source flow, W* are the path flow-rates, Ry i = 

are the valve resistances, and 6i are the angular valve positions. 

2) Actuator Submodel: Each flow path has an actuator 
valve unit consisting of a ball valve, brushless dc motor, and 
built-in servo speed controller. The angular valve speed (w*) 
response of this unit to an input valve speed command (u{) 
is given by the actuator vendor [5] as a first-order lag. This 
is represented as 



The resulting valve speed is integrated to achieve the valve’s 
position, 'fhis is modeled as 

Oi = Wi, Vi = M- (3) 


w* = - ( k v Ui - Wi), 
T 


Vi = 1,4. 


The simple actuator model represented by (2) and (3) is an 
acceptable model for this problem because the actuator is not 
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Fig. 3. TCS tracking profiles. 


explicitly controlled, but only used as an effector whose gain 
and phase responses are adequately modeled. 

3) Pump Submodel : The centrifugal pump model [5] con- 
sists of a pump motor, pump speed controller, and system 
pressure controller equations 


robust, and decoupled tracking performance of time-varying 
flow-rate profiles W* (t) is desired. The problem is defined 
mathematically as 

lim ei(£) = 0, Vz = l,3 (5) 

t—*OQ 


u p = - — ( k t I a - kiJ 1 ) 

3 ml 

f a — [ Rrfila k e UJp k 2 P e H~ QX C ] 

- — ( w p - w /) 

T f 

X c =MX c + B-\u f P P‘] T (4) 

where u p is the pump speed, I a is the pump motor armature 
current, utf is the pump speed feedback into the pump speed 
controller, X c is the 2 x 1 vector of system pressure controller 
and pump speed controller states, and P e is the system pressure 
error (P* - P). 

4) TCS Model Summary: The composite 19th-order multi- 
input-multioutput nonlinear TCS with parameter uncertainties, 
interaction disturbances (dynamic coupling between parallel 
paths), and external disturbances is represented by (l)-{4). The 
model contains three input valve commands (ui,u 2 , and u$) 
and three output flow-rates {Wi,W 2 , and W$). The model 
is verified, through comparisons with test data from Space 
Station TCS tests, to adequately represent the salient dynamics 
of the SSFF TCS. 

B. Control System Objective 

The control system objective is to simultaneously control 
the thermal environments of several furnaces and related sub- 
systems plumbed in parallel. Specifically, the high accuracy, 


where 


e i (t) = W*(t)-W i (t). (6) 

The solution of this problem using sliding mode is two fold 
[6]. First, sliding surface equations 

<7i(x) = 0. Vi = 1,3 (7) 


are chosen such that the flow-rate tracking errors defined by (6) 
are the outputs of linear homogeneous systems with constant 
coefficients and with desired eigenvalues placements. Second, 
discontinuous control laws of the form 


„.(+ _ / uf{t,x) (Ti(x) > 0 

1 X (t,x) (Ti(x) > 0 


Vi = 1,3 (8) 


are chosen such that the desired output tracking motion goes 
to and stays on the sliding surfaces defined by (7). 


III. Sliding Mode Controller Design 

This section discusses the theory of sliding mode control as 
it is applied to the design of the TCS problem. This section 
develops the system transformation that is most useful in de- 
termining the system’s vector relative degree and in analyzing 
the stability of the system’s internal dynamics; synthesizes the 
sliding surfaces upon which the desired output responses are 
restricted; and develops the sliding mode controllers that force 
the system to move and to stay on the sliding surface for all 
subsequent time. 
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Fig. 4. Normalized external disturbance. 



Fig. 5. Path 1 flow tracking response. 


For the TCS output tracking problem, the linear homoge- 
neous sliding surface equations are defined as 

r, -1 

= ^2 c ji e i 3 \ V? = 173 (9) 

3=0 


where the relative degree r* for each output channel and the 
system’s internal dynamics, if any, are determined through a 
system transformation. It is well known [7] that sliding mode 
control is applicable in systems with stable internal (zero or 
asymptotic zero) dynamics. 
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Fig. 6. Path 2 flow tracking response. 



Fig. 7. Path 3 flow tracking response. 


A. System Transformation 

The 19th-order (n = 19) TCS system represented by (1) — (4) 
is presented in matrix form 

x = f(x) + Gu 

y=h{x) (10) 

where x is a 19 x 1 vector of model states, f(x) is a smooth 
19x1 vector field, G is a constant 19x3 matrix, u is a 3 x 


1 vector of input commands, y is a 3 x 1 vector of outputs, 
and h(x) is a smooth 3 x 1 vector field. 

The system transformation begins by deriving a direct 
input-to-output relationship. This is accomplished by taking 
successive derivatives of each output until at least one input 
appears [7] 

m 

2 /< r,) =L r jh t (x) + J2 LgL r j~ 1 h t (x)u j , = m 

j=i 

(ID 
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Fig. 8. System pressure response. 


where the scalar function Lrh is the Lie derivative [7] of 
h with respect to /, m is the number of inputs, r* is the 
number of differentiations required for the input the appear 
(called the ‘‘relative degree” of the zth output channel), and 
Lg.L T j~ l hi{x)uj 7 ^ 0 for at least one j. Applying ( 11 ) to 
(10) for m = 3 yields the ninth-order (r = 9) input-output 
description 


V l 
y 2 
y' s 



rL Pl 

a0W 2 e~^uj 2 

tLp 2 

a/3W*e-e e *uj 3 
tLp 3 


+ D(x) 


u i 
u 2 

^3 


( 12 ) 


where we have (13), shown at the bottom of the page, and 


V u = -a0 2 u 2 i W?e- fie ' + 


2Q li W? + 2n ll W l 

L Pi 


a/3u>ie 001 

(14) 


m _ ! 2v0WiWiSlueP ei ( Oit ( A TJy r2r»2 , OA n \ 

^2t — T 4* 2 (41 4" 2f2iif22t) 


*3i " 


Lp, 

2 n u 
T~c 

4 


L Pi C 


n 2l w t w s -J 2 Wj 


j = i 


*«=*• 




k^x 


«ii ----n 2t w? + p 

ft 2i ® = -{ae~ 0ei + Ki), Vi = 175. 


(15) 

06) 

(17) 

(18) 
(19) 


Matrix i)(x) is full rank and nonsingular for any operational 
mode of tie system. The system (12) has a vector relative 
degree of - = [3,3,3] and a total relative degree of r = 9. 
Since r < (n - r)th or lOth-degree internal dynamics exist 
and shoulc be analyzed for stability. 

The system (12) is then transformed to the phase- variable 
canonical torm [ 8 ] and combined with the internal dynamics 
[dynamics not involved in ( 12 )] to form the so-called “normal 


D(x) 


aPKvW? 00j 

rL Pl 

0 


0 

o0KyWl _ ge , 
rh P2 6 


0 

tL P3 


(13) 


0 


0 
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form’' representation of the original system 

Z\ x — Z2i 
Z 21 — 

1 ^ 1 

Z 3t =— Y 9ji + - af3W?e ~ 0e ' (<* + K v m) 

Lp - [U 

Vi = Zu, V* = 1,3 
Vk =qk{Z*T)), V/c = 10.19 (20) 

where W x . 6 x . and uj z can be expressed through Zu,Z<zi, 
and Z 3t using (14)— (19). The internal dynamics 77 consist of 
equations in (1H4) that represent the system pressure drop, 
pump source flow, pump, and path 4 submodels. 

Stability of the internal (asymptotic zero) dynamics is shown 
using Lyapunov’s first method [10]. This method involves 
linearizing the nonlinear internal dynamics throughout the 
operating range of the system. Eigenvalues at each operating 
point are calculating to assure that all have negative real parts 

mi 

B. Sliding Surface Design 

The second major step in the design process is to choose the 
desired linear homogeneous sliding surfaces upon which the 
motion of the tracking errors ( 6 ) is restricted. According to 
(9), where r x = 3 ,Vi = 173, the sliding surface equations 
are chosen as the set of second-order linear homogeneous 
differential equations with constant coefficients of the form 

<Ti(x) — e‘i -b cuCi + c\iCi = 0, Vi = 1,3. (21) 

The coefficients are chosen to obtain the desired transient 
response of the tracking error motion. These coefficients are 
chosen as cu = 225.0 and c<n = 21.0, Vz = 1,3, using the 
ITAE method [12]. According to the standard second-order 
form 

e + 2( ) uj n e + — 0 (22) 

the sliding surface has a damping ratio of £ = 0.7, a natural 
frequency of = 15.0, and a settling time of 4 /£u; n — 
0.38 s. This represents a desired transient response of the 
system operating in sliding mode; however, there are several 
realization issues (for example, discrete sampling rate and 
unmodeled modes) that must be considered before this desired 
transient performance is obtained. Section IV will address 
these issues. 

C. Control Function Design 

The third major step is to design the control function to 
provide the existence of sliding mode on the sliding surfaces 
(21). The motion of the system in a-space is described by the 
system 

*£>[/(,) + <** (23) 

A candidate Lyapunov function [9] is introduced as 

V{x)=\ cj T (x)a(x) 


where V(x) is positive definite and 

lim V (x) = 00 . (25) 

IK*)IHoo 

Global asymptotic stability of a(x) = 0 is assured if the 
derivative of (24) is negative definite [9]: 

V(x) = cr T (x)a(x) < 0. (26) 

Moreover, if the inequalities 

<Ti(x)&i(x) < -pi\oi{x)l Pi> 0. Vi = 1,3 (27) 

are satisfied, then finite sliding surface reaching times are 
calculated as 

^ -JfiM \/i = 173. (28) 

Pi 

Analysis of the system (23) shows that the matrix 
[da(x)/dx\G is diagonal. This allows system (23) to be 
rewritten in decoupled form 

di(x) = di(x) - bi(x)ui, Vi = 1,3 (29) 

where 


bi(x) = 


+ C2i&i + Cue 

aBk v W?e~ 08 ' 


apUiWfe-W' 


(30) 

(31) 


and bi(x) > 0. Vx. The discontinuous control law 

Ui(x') — Ui sign[(7j(x)], V* = l,3 (32) 


where 

Ui>u\\x)+Ji- (33) 

bi(x) 

«r(x)=[6,(x)]" 1 di(x) (34) 

satisfies inequality (27), Vi = 1,3. This means that the control 
law (32) provides the existence of sliding mode on the sliding 
surfaces (21) with finite reaching times (28). 

The equivalent control (34) is calculated through a typical 
flight simulation, including disturbances, and shown to be 
mostly contained within a ±5 volt boundary. Considering this 
boundary and that the actuator is designed to accommodate 
input control signals between ±5 V, Hi is chosen as 5.0 V, 
and the resulting discontinuous control function (32) is shown 
as . 


(24) 


Ui{x) = 5 - sign[cri(x)], Vi = 1,3. 


(35) 
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Fig. 9. Discrete sliding control input to path 1 . 



Fig. 10. Sliding surface response in path 1. 


IV. Realization of Sliding Mode 

The SSFF computer samples measurements once every 
T = 0.005 s and utilizes 12-bit A/D converters to digitize 
the data. It can be shown that the discrete sampling time 
T = 0.005 s satisfies the sufficient condition for the existence 
of sliding mode for systems in discrete time [14] 

|a < [*r]|<|a < [(fc-l)T ]| 1 = 


The digit! *ed signals are low-pass filtered using a fourth-order 
digital butterworth filter (with a passband of w c = 100.0 
rads/s) to remove measurement noise [13]. The flow tracking 
errors ei[ zT ] are then computed and combined to form the 
sliding surfaces (21) in discrete form 

a t [kT - c\ei[kT] + c 2 ei[(k - 1 )T\ + c 3 ei[{k - 2 )T], 

Vi = 173 (37) 


(36) 




620 


IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 6, NO. 5, SEPTEMBER 1998 



Fig. 11. Source pump map. 


where k = 


0, 1, 2, ■ • • , and 

225. OT 2 + 21. 0T + 1.0 


Cl 


c 2 = 


T 2 

21.0T+ 2.0 


1 

c 3 = Y 


(38) 

(39) 

(40) 


In general, sliding performance (specifically bandwidth) can 
be limited according to the discrete sample time and any 
unmodeled delays in the system [7], This is shown as 


u) n < min 


27 r 27T 

5 T 1 3T4 


(41) 


where u; n is the sliding surface bandwidth, T is the discrete 
sample period, and Ta is the largest unmodeled delay. For the 
TCS model, the two largest unmodeled delays are through the 
actuator (T^ l < r) and the noise filter (Ta 2 = 27r/u; c ). The 
delay Ta 2 through the noise filter is calculated as the dominate 
delay. Solving inequality (41), it is shown that the prescribed 
/ ' e bandwidth (oj n = 15) specified in (37) is not 

iinuted and is thereby realizable. 

key issue in the realization of sliding mode is 
xY wCrH]-.:*iOn of chattering — defined as the phenomenon 
of nonidc high speed switching. Defining the “smoothing” 
sliding mode controller 


Ui(kT) 

( 5 -sign[ai(kT,x)], \<n{kT,x)\ > e 
~ | - cri(kT,x), |cr i(kT,x)\<e, Vi = 1,3 

( 42 ) 


smooths the discontinuous control function (35) within e- 
vicinity of the sliding surface, thereby eliminating chattering. 
The boundary width e = 1.5 is chosen to minimize steady-state 
error (integral control may be added if this is significant) and to 
achieve a tradeoff between system robustness and chattering. 
This tradeoff occurs because ideal nonlinear sliding control 
is robust (or invariant) to disturbances and uncertainties, 
whereas the addition of a linear region (smoothing effect) 
causes the system to lose some invariance to disturbances and 
uncertainties within the specified boundary width. 


V. Simulations 

The mathematical parameters for the numerical simula- 
tion of the TCS model are defined in the Appendix. The 
Runge— Kutta fourth-order algorithm with a 0.001 -s step size 
is used for the integration of the plant differential equations, 
and the discrete sliding mode controllers are scheduled every 
0.005 s. The performance of the sliding mode controllers is 
assessed under three operating conditions: 1) simultaneous 
flow profile changes in each of the three parallel paths; 2) 
external disturbance torques applied to the flow actuators; 
and 3) parameter uncertainties occurring in the model. All 
measurements include uniformly distributed white noise of 
amplitude 0 to 1 lbm/hr. 

Typical flow tracking profiles for paths 1, 2, and 3 are 
defined in Fig. 3. In response to these profiles, each path’s 
sliding mode controller modulates its flow control valve to 
track the profiles. Parallel valve modulations cause the system 
pressure to deviate from its set-point P* 9 therefore causing 
interaction disturbances to each path. The pressure control (a 
decoupling effect) is considered poorly designed or operat- 
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in g in a degraded mode to intensify the interaction distur- 
bances. 

External disturbances such as unwanted torques to the 
control actuators originate from a variety of sources such 
as power delivery spikes, electromagnetic interference, and 
flow turbulence variations. Such torques are unmeasured and 
can greatly affect the performance of the TCS. For sim- 
ulation purposes, torques are normalized (see Fig. 4) and 
added to the control input signal to represent an external 
disturbance. 

Plant uncertainties are considered because of the difficulty in 
accurately modeling the system under investigation. Uncertain- 
ties can be the result of imprecise measurement instruments, 
improper assumptions, incorrect measurement procedures, or 
other measuring constraints. For this simulation, the source 
resistance Rs , the path flow inertance Lp x , and the valve gain 
k V2 are considered 120, 80, and 80%, respectively, off their 
nominal values to represent possible plant uncertainties. 

The resulting flow-tracking responses in paths 1, 2, and 3 
are shown in Figs. 5-7, respectively. The solid and dashed 
lines represent the responses to the sliding mode and PID 
controllers, respectively, and the dotted lines represent the set- 
point profiles. Seen in Fig. 8, the system pressure deviates 
from its set-point significantly for both the sliding mode and 
PID simulations. This deviation can be a disastrous effect 
to controlling flow in the parallel paths. Nevertheless, the 
sliding mode response tracks the desired profile almost exactly, 
whereas the PID response deviates from the desired profile 
by as much a 150%. This loss of performance would result 
in the SSFF having to shut-down; however, implementation 
of sliding controllers would increase system reliability and 
maintain system performance. 

The sliding control input (ui) shown in Fig. 9 is smooth 
(does not exhibit chattering). The sliding motion (< 7 i) shown 
in Fig. 10 is within the 6-vicinity of the sliding surface except 
for a few short intervals, where lack of control resources causes 
sliding mode to be destroyed. In summary, robust precise 
decoupled tracking performance is demonstrated throughout 
the operating range of the TCS, even in the event of system 
pressure control failures. 


APPENDIX 

Mathematical Model Parameters 


Hydrau ic Submodel Data. 

Path Resistance, Rp = 3.57 x 10~ 5 Wi [lb / hr 2 ] [lb m 
> 2 ]- x - 

Path Inertance, L Pi = 7.03 x 10 4 [hr ■ s • lb/] [lb m in 2 ]* 1 . 
Valve Loss Coefficients: 

a = 0.129 [lb / hr 2 ][lb^ in 2 ]- 1 
(3 — 0.169 deg" 1 . 

Source Resistance, Rs = 1.0 x 10~ 5 [lb / hr 2 ][lb m in 2 ]- 1 . 
Source Inertance, L$ = 1.1 x 10~ 3 [hr ■ s * lb/][lb m hr] -1 . 
System Compliance, C = 0.1 [lb/ hr][lb m nr s]" 1 . 

Source Pressure, Ps = f(ui p: Ws) See Fig. 11 for map. 

Actuator Submodel Data: 

Valve Time Lag, r = 0.01 s. 

Valve Constant, k v = 1.0 [deg] [s * v]" 1 . 

Pump Submodel Data: 

Rotor plus Load Inertia, j m i = 3.1 x 10 -5 in • lb/ s 2 . 
Motor Gain, k t = 0.33 [in • lb/] [amp]* 1 . 

Motor Inductance, Lm = 0.0012 henrys. 

Motor Resistance, Rm = 0.618 Q. 

Speed to Torque Gain, k\ = 5.07 x 10“ 7 [lb / s 2 in] [rad 2 ]* 1 . 
System Pressure Set-point, P* = 10.0 psi. 

Pump Speed Controller Feedforward Gains: 

&2 = 0.0021 [amp] [rad] -1 . 
fc 3 = 0.021 [amp- in 2 ] [lb/ s]" 1 . 

Pump Speed Lag, r/ = 0.008 s. 

Motor Lack EMF Gain, k e = 0.042 [v][rad • s] _1 . 

System Pressure and Pump Speed Controller Matrices: 


Q = 0.0213 400.0] 


M — 


0 

4.3 x lO" 4 


0 

0 




0 

-4.3 x 10“ 4 


- 0.1 - 0.1 
-4.3 x 10* 5 4.3 x 10~ 5 
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VI. Conclusions 

Sliding mode control theory has been employed in the 
context of the TCS problem. The discrete realization of the 
controllers has been developed. The most significant conclu- 
sion from this work is the determination that sliding mode 
control provides precise, robust, and decoupled control of flow 
profile tracking in the presence of interaction and external 
disturbances, parameter variations, and system nonlinearities. 
Sliding mode provides this control over the full operating 
range of the TCS; moreover, this control is maintain when the 
system pressure control is operating in a degraded mode. Ex- 
cellent simulation performance and robustness is demonstrated 
for the discrete sliding mode controllers. Implementation of 
the sliding mode controllers will enhance the SSFF system’s 
reliability and performance without additional cost or resource 
impacts to the system. 
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